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\^ '. Abstract 

^ ■ We present the first direct determination of the running top-quark mass based on the total cross sec- 
O ■ tion of top-quark pair-production as measured at the Tevatron. Our theory prediction for the cross 
section includes various next-to-next-to-leading order QCD contributions, in particular all loga- 
i rithmically enhanced terms near threshold, the Coulomb corrections at two loops and all explicitly 
5-j I scale dependent terms at NNLO accuracy. The result allows for an exact and independent variation 
of the renormalization and factorization scales. For Tevatron and LHC we study its dependence 
on all scales, on the parton luminosity and on the top-quark mass using both the conventional pole 
mass definition as well as the running mass in the MS scheme. We extract for the top-quark an MS 
mass of m(ju = m) = 160.0^2 2 ^s^' which corresponds to a pole mass of = 168.9^2 4GeV. 



The top-quark is the heaviest elementary particle discovered so far and it is likely to be the most 
sensitive probe of the electroweak symmetry breaking. This is reflected in the fact that in many 
extensions of the Standard Model the top-quark plays a special role. The precise measurement 
of top-quark properties is thus an important task for the Large Hadron Collider (LHC) (see e.g. 
Refs. [1, 2]). One of the most basic quantities in that respect is the total cross section which is 
currently measured at the Tevatron and will be measured at the LHC. The precise measurements 
aimed for at the LHC are asking for an equally precise theoretical prediction to compare with. 
In this Letter, we update and extend the predictions of Refs. [3,4] for the total hadronic cross 
section of top-quark pairs and its associated theoretical uncertainty. Related recent studies have 
also appeared in Refs. [5-7]. As a novel aspect of this Letter, we employ the MS definition for the 
top-quark mass and present the total cross section as a function of the running mass. This allows 
the direct determination of an MS mass from Tevatron measurements for the total cross section [8]. 

We start by recalling the relevant formulae for the total cross section of top-quark hadro- 
production within perturbative Quantum Chromodynamics (QCD): 

J s 

o-pp^ttx(S,mf) = ^ 2] I ds Lij{s,S,nj) fij{p,M,R), (1) 

S ^ 

2 \ r ds I ^ 2\ ( ^ iX 

Lij(s,S,fif) = - I -j'^i/p\j^l^f\<Pj/py^^l^fj^ (2) 

s 

where S and nit denote the hadronic center-of-mass energy squared and the top-quark mass (here 
taken to be the pole mass), while Ljj is the usual definition of the parton luminosity with the 
parton distributions (PDFs) ^i/p at the factorization scale /i^. The scaling functions fij parame- 
terize the hard partonic scattering process. They depend only on dimensionless ratios of mt, fip 
the renormalization scale yti^ and the partonic center-of-mass energy squared s, with the defini- 
tions p = 4mj/s, R = /J-r /fij- and M = fij/rn^. The radiative corrections to the scaling functions 
fij at the next- to-leading order (NLO) [9-1 1] are long known and recently even the complete ana- 
lytic expressions have become available [12]. In order to quantify the theory uncertainty also the 
next-to-next-to-leading order (NNLO) must be included. Presently, these NNLO corrections [3,4] 
are approximated by the complete tower of the Sudakov logarithms and including all two-loop 
Coulomb corrections. 

In the present study we include consistently the channels, qq, gg and gq through NNLO and 
we provide parametrizations for all necessary scaling functions in the standard MS scheme for 
mass factorization. This allows for an easy handling in phenomenological applications. Our 
phenomenological study reflects the latest measured value for the top-quark mass [13], nit = 
173.1;^} 3 GeV, and employs new PDF sets [7, 14]. Let us briefly summarize the key aspects of 
our update with respect to Refs. [3,4]: 

i. ) We use exact dependence on the renormalization and factorization scale. This allows for 

an independent variation of yu^ and //^ (extending Ref. [15]) and is commonly considered 
as more reliable to establish the theoretical uncertainty of perturbative predictions (see e.g. 
Ref. [16]). 

ii. ) We perform the singlet-octet color decomposition consistently when matching our threshold 

expansion at NLO using results of Refs. [17-19]. The numerical impact on phenomenology 
at LHC and Tevatron turns out to be negligible, though. 
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iii. ) We discuss those residual systematical uncertainties of our predictions for LHC and Teva- 

tron, which are inherent in the approach based on threshold resummation and we comment 
on the size of unknown corrections. 

iv. ) We quantify the numerical impact of other known effects on the total cross section, such as 

the electro-weak radiative corrections at NLO [20-22] and bound state corrections in QCD 
at threshold [18,19]. 

V.) We study the dependence of the total cross section on the definition of the mass parameter. 
For the conversion of the conventionally used pole mass nif (see Eq. ([T])) and the scale depen- 
dent MS mass mip.^) we exploit well-known relations to NNLO [23] (see also Refs. [24,25]). 
We investigate the apparent convergence of both definitions, nit and mip.,), in perturbation 
theory through NNLO. 

We also employ the analytic results for the NLO scaling functions [12]. As a net effect these lead 
to small improvements in the gq- and the ^g-channel contributions of our NNLO prediction. 



Theoretical setup 

The perturbative expansion of the scaling functions fij in the strong coupling Us up to two loops 
around M = R = 1, i.e. mt= reads 



(3) 



and the functions fy , fy and fy contain, at each order in a^, genuinely new information to be 
calculated from first principle in perturbation theory. At Born level. 
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where J3 is the heavy quark velocity with = ^Jl -p. At NLO the known functions fj^j^^ can be 
described through parametrizations which are accurate at the per mille level. Our one-loop fits use 
the following ansatz. 
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where n f denotes the number of light quarks and we have kept the complete dependence on in all 
parametrizations manifest. The Sudakov logarithms \n/3 at threshold and the Coulomb corrections 
(~ 1 //3) in Eqs. (IT])-® are exact [9]. The constants read 

^ i?^_JL,_Hl!!il,15!£!2. 0.03294734, (10) 



324;r 1296 108 ;r 21 n 
nil 283 89 ln2 7 ln22 
230471- ~ 18432^" 128 n ^T6~ 



= ^z:z7r. 777:^^-^ + ^ « 0.01875287. (11) 



They are known from the computation of the NLO QCD corrections to hadro-production of quarko- 
nium [17] (see also Refs. [18, 19]), where also details of the decomposition of a^^ for color- 
singlet and color-octet states can be found. The constants a'^ in Eqs. ([TO]) and ([TT]) emerge from 
Refs. [17-19] by means of a simple Mellin transformation and agree with the values quoted in 
Ref. [12]. The coefficients of the functions h(J3,ai,...,ai'j) and h^gqifi,a\,...,ais) in Eqs. dT])-® 
are determined in a fit to the analytic expressions ofRef. [12]. Near threshold we have /7(/3) =0{fi^). 
More details are given in Eqs. (IA.7I) and (|A.8I) . Due to the larger number of parameters in the fit 
functions, it is evident that Eqs. (|7])-® supersede earlier parametrizations [9] with respect to ac- 
curacy. 

At two-loop level we know the complete tower of Sudakov logarithms, In^jS with k= 1 , 4, for 
the functions f*^^^ and f'^g^ and, in addition, also the complete Coulomb contributions, ~ 1 /jS^, 1 //?. 
The channel gq is power suppressed near threshold relative to qq and gg. However, extending soft 
gluon resummation to power suppressed quantities (see e.g. Refs. [26, 27]) and using Eq. ([8]), we 
can determine (at least) the leading term ~ Xvc' [5 of the function f'^q''- We find 
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Eqs. ([T2]) and ([141) are exact up to the unknown constant terms C^?* and C^^^ of order O{0^), whereas 



Eq. ([T3]) receives further corrections of order Oifi^ In^yS). Please note that the numerical coefficients 
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in Eqs. (fT2l) and (fT4l) have slightly changed compared to Ref. [3]. The ln^/3-terms in Eqs. (fT2l) and 
(fT4l) are affected by using the exact coefficients (flOl) and (fTTI) in the matching at NLO. The linear 
terms proportional to In^S in Eqs. (fT2l) and (fT4l) contain genuine two-loop contributions. Among 
those is the soft anomalous dimension D^qq (see Ref. [3]). Inserting the respective numerical coef- 
ficients we find 730.73916 + 23.776275n/ in Eq. ^ and 3035.5764-0.88761378n/ in Eq. ([T?]). 
The latter value has changed with respect to Ref. [3] due to a consistent separation of the color- 
singlet and color-octet contributions in a^^ at NLO [17-19]. However the phenomenology is rather 
insensitive to this change and is only aff'ected at the per mille level. The Coulomb terms (~ 1 IP) in 
Eqs. (fT2l) and (fT4l) contain all contributions from the two-loop virtual corrections. Eq. (fT3] ) gives 
the leading (though formally power suppressed) contribution at two loops to the ^(^r-channel. We 
include fg^^ in our analysis for three reasons. Firstly, under evolution of the factorization scale 
the g^-channel mixes with the two other channels and for a consistent study of the factorization 
scale dependence this channel also needs to be taken into account. Next, the luminosity Lgq in 
particular at LHC is sizeable and Eq. (fT3] ) ofl^ers a way to control its numerical impact at higher 
orders. Finally, Eq. (fT3] ) provides a first step towards a general study of power suppressed but 
logarithmically enhanced terms near threshold for top-quark production. 

In Eq. ([T]) the dependence of the scaling functions fij on the renormalization and factorization 
scales, and ix^ can also be made explicit. Starting from the expansion in through NNLO 
around R=l, i.e. yu,. = yu^, we introduce 

fijipMA) = flf{p) + Anas[fl]''^{p) + LMfl]'^{p)] 

+(4;ra,)'|i;.f V) + ^M/;f V) + 4i:f ^(P)) ' (15) 

where we abbreviate Lm = ln(nj./mj). The logarithmic tower in M = fxj./mj, that is all terms 
proportional to Lm in Eq. (fT5l) can be derived by renormalization group methods in a straight 
forward manner. The explicit results in the MS scheme for all channels read at NLO 

= -(2P^I^^-2/3om4j\ (16) 

.(11) _ _p(0)^.(0)_^p(0) .(0) .,7. 

where ij = {qq,gg} in Eq. (fT6l) . At NNLO we have 

41" = -(2n';-2/*.l>)«4°'-a<'-3All)«C-2''i>/r. (18) 

4r = ~<'®c-<n>fi?-3/soi)«/r-4'>4°'-2^n^4?' (20) 

A? = ^/■™«(3<,' + ''S'-5Al)«/,T4P^>(n? + 34?-5An)«4?. (21) 

fS" = -(2P^^2/)a)®/™-(2^■s'-3All)»/^-2n>/^. (22) 

4f = 5i^<><>4°U(<><^24>pg'-50„ + 3^^1)«4°'. (23) 

In Eqs. (fT6l)-(l23]) the (S) products have to be understood as standard convolutions and a sum over 
all active quarks and anti-quarks is implied as well. The coefficients of the QCD jS-function are 
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The splitting functions P. . can be taken e.g. from Ref. [28, 29]. At leading order they read 



(167rVSW = -( — -2-2x + 35(l-.^)), (25) 
(16;rVgg(-^) = 2n/(l-2x + 2jc2), (26) 
(16;rVS^(x) = l(l-4 + 2x), (27) 



3\x 



\ -1- XX J / ^ 



where the factor {I6n^) accounts for the normalization used in Eqs. (fT6l)-(l23l). In general, we 
have Pfjithis article) = (\6n^y^'^^^ P\'j(Refs. [28,29]). Please also note the explicit factor of (2n/)"^ 

in Eqs. (fT7])- ((23l) . which is due to the definition of Pf^ in Eq. ^ and P^j,^ in Ref. [29]. Simple 
fully analytic expressions for fq^^\ fgq^^ and fgg^^ are long known [9] and precise fits for all 

scaling functions fij^\f^ in Eqs. ([T8])- (|23]) . typically to per mille accuracy, are presented in 
the Appendix in Eqs. (|A.1I) - (|A.6I) and Tabs. \5\^ Finally, the complete scale dependence for 
fij{p,M,R) in Eq. ([H) with ix^. 4^ ix^ is easily obtained as 

f,jip,M,R) = Aj(p,MA) + 4nas[2l3oLRf^^} (29) 

fgq{p,M,R) = fgq{pMA) + {4nasf[3liQLRfl'^K3liQLMflq^], (30) 
where Lr = \n(pi^/iuj) and ij = {qq,gg} in Eq. (|29l) . 



Phenomenological applications 

We are now in the position to address the phenomenological consequences. The approximate 
NNLO prediction which includes exact dependence on all scales is based on Eqs. (IT2l)-(IT4l). (fT8l) - 
(l23l) and (|29l ). (I30l) . If not otherwise stated, the top-quark mass is the pole mass at = 173 GeV. 

Let us start with the scale dependence of the NNLO cross section as shown in Fig. [H Our 
study of the theoretical uncertainty allows us to assess the effect of independent variations of the 
renormalization and factorization scale ju^ and //^ in the scaling functions for the hard partonic 
scattering process in Eq. ([T]). In doing so, one should keep in mind however, that all currently 
available PDF sets from global fits always fix /i^ = /i^. We define the theory uncertainty arising 
from the independent variation of ju^ and /i^ in the standard range fir,f^^ 6 [mt/2,2mt] as 

min cr(ju^,fijr) < cr{mt) < max cr{n,.,Hj:). (31) 

The contour lines of the total cross section for the LHC and Tevatron arising from this procedure 
are shown in Fig.[T] The standard range /i,.,//^ e [mf/2,2mf] corresponds to the region displayed in 
Fig.[T]because of logjo(2) ~ 0.3. We have normalized all results to the value at yu,. = //^ = and the 
variation with fixed scales //,. = /i . proceeds along the diagonal from the lower left to upper right in 
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Figure 1 : The contour lines of the total hadronic cross section from the independent variation of renormal- 
ization and factorization scale //^ and jij, for LHC with = 14TeV (left) and Tevatron with yfs = 1.96TeV 
(right) with MSTW 2008 [14]. The cross sections are normalized to the values at yu,, = fi^. = nit and the range 
corresponds to yu^,jUy € [mt/2,2m,]. 



the plots. For Tevatron in Fig. |l(b)| we see, that the gradient in the complete (/i^,//^)-plane is almost 
parallel to this diagonal, thus the uncertainty according to Eq. (|3TI) remains -5% < Act < +3% in 



very good agreement with previous results [3,4]. At LHC in Fig. |l(a)| the maximal deviations in the 
Ou^,/Zy)-plane are -3% located at (2m,, 2m/) and +1% at roughly {mt/2,2mt), so that the uncertainty 
range (|3TI) becomes slightly larger, -3% < Act < +1% compared to what has been derived before 
with a fixed scale ratio ju^ = jUjr. Very similar numbers for both colliders are obtained with the 
PDF set CTEQ6.6 [7]. Also recall, that we include the gt^-channel through two loops. Thanks to 
Eqs. (fT8])-(l23]) we control the exact scale dependence also for all contributions proportional to the 
parton luminosity Lgq. We conclude from Fig. [T] that the theoretical uncertainty due to ju^ and 
variation is well estimated by case of identical scales /u^= /u 



f 

In order to quantify the PDF uncertainty we apply the standard definition 



^ 1 E (o-t^-o-i,-)^ 

\^k=l,npDF 



(32) 



which determines Act from the variations crk± of the cross section with respect to the k parameters 
of the PDF fit. Typically the PDF error is added linearly to the theory uncertainty obtained from the 
scale variation. This is the commonly adopted choice and we employ the PDF sets CTEQ6.6 [7] 
and IVISTW 2008 [14]. The latter set gives two uncertainties at different confidence levels (CL), 
one at 68% CL and the second at 90% CL. Throughout this study we use 68% CL only. IMoreover 
in the chosen interval fi e [mt/2,2mt] for a given fi = ju^ = /u^ the error ActQj.) in Eq. (|32|) has only 
a very weak scale dependence. That is to say we find to good accuracy AcrQi = nit 12) - AcrQi = 
ntt) - A(t(iu = 2mt) so that a determination of Act at the central scale /u = nit should suffice for all 
practical purposes. 
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Figure 2: The mass dependence of the total cross section at NLO (green) and approximate NNLO (blue) or- 
der for LHC at VS = 14TeV (left) and Tevatron at - 1.96TeV (right) and the PDF set MSTW 2008 [14]. 
The bands denote the theoretical uncertainty from scale variation keeping //^ = ju^. and the PDF uncertainty 
in the range [mt/2,2mt]. In Fig. |2(b)| for the Tevatron the NLO and approximate NNLO bands overlap only 
partially giving rise to the corresponding medium band (dark green). 



In Fig.[2lwe show the mass dependence of the total cross section, comparing the NLO and our 
approximate NNLO prediction. The band summarizes the total theoretical uncertainty from the lin- 
ear combination of the scale uncertainty for the case //^ = //^^ and the PDF uncertainty Eq. (|32|) . We 
display the LHC and Tevatron predictions using the MSTW 2008 PDF set [14]. The improvement 
of the NNLO prediction is manifest for both colliders. 

Next, we discuss the sources of remaining systematical uncertainties. Undoubtly, a complete 
calculation of the complete NNLO QCD corrections to hadronic top-quark pair-production would 
be highly desirable (see Refs. [30-37] for progress in this direction). This lacking, the main sys- 
tematic uncertainty of our approximate NNLO result are the sub-leading terms in the scaling func- 
tion fl'^'^\ They might become accessible by extending the approach of Refs. [26, 27] or else 
could ^e modeled through power suppressed terms in Mellin space (see e.g. the scheme A = 2 in 
Ref. [38]). By including the leading term for fg^^^ we have taken a first step in this direction and 
we have found numerically small effects only. In order to quantify our systematical uncertainty, 
we adopt the following prescription: We compute the ratio (TNLL+Coui/cexact at one loop, where 
o"NLL+Coui only contains the Sudakov logarithms and the Coulomb terms in /^i^^ and fgl^\ i.e. 
the content of the square brackets in Eqs. (|7]) and dH). This checks how well the exact hadronic 
cross section in Eq. ([T]) is approximated, if only the threshold approximation enters in the con- 
volution with the parton luminosities. Typically we find crNLL+Coui/o"exact ~ 0.7(0.9) for the LHC 
(Tevatron). If translated to the genuine two-loop contribution (see e.g. Tab. [T] below), then a sys- 
tematic uncertainty of 0(30%) implies a cross section uncertainty of Act - (9(15) pb at LHC and 
of Act - 0(0.2) pb at Tevatron. These numbers are corroborated by other observations, like the 
generally small impact of the ^(^-channel which is entirely sub-leading. 

How does this affect the previous discussion of the scale dependence? Let us define a system- 
atic uncertainty Asys obtained from a variation of the scaling functions fj^J^^ in Eqs. ([T2l)-([T4l) by 
+30%. All other scaling functions are known exactly anyway. The result for our NNLO cross 
section (always normalized to the value at ju = nij and Asys = 0) is shown in Fig. [3]for // = /i^ = /i , 
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Figure 3: The scale dependence of the approximate NNLO cross section ctnnlo for the choice ji-n,.- fi^. 

using the MSTW 2008 PDF set [14] at LHC with y/S - 14TeV (left) and Tevatron with Vs = 1.96TeV 
(right). Asys denotes the estimated systematic uncertainty of our threshold approximation at NNLO and all 
results are normalized to the value of cr(ju = nit). 



and the MSTW 2008 PDF set [14]. It is obvious that the predictions are very stable within the 
standard range yu 6 [m//2,2m;] for all cases, i.e. A^ys = and (T + Agyg. For the case A^yg — we find 
a variation of -3% < Act < +0.5% for LHC and -5% < Act < +3% for Tevatron (compatible with 
Fig. [I]) and similar numbers for the other two cases, cr + Agys. 

In the present analysis, we have also neglected the effect of the new parton channels qq, qq 
and qtqj (for unlike flavors / ^ j), which come in through real emission at NNLO only. Important 
insight can be gained here from the recent calculation [39,40] of the radiative corrections for tt+l- 
jet production at NLO, because they represent a significant part of the complete NNLO corrections 
for inclusive top-quark pair-production. At NLO tt+l-jet production contains the one-loop one- 
parton real emissions as well as the double real emission processes, and the latter also include the 
above mentioned new channels. It was found that the radiative corrections at the scale /i,. =// . = 
are rather small. Depending on the kinematical cuts (e.g. on the transverse momentum of the jet) 
they amount to (9(20) pb at LHC and to 0(0.2) pb at Tevatron (see Ref. [40]). This provides further 
evidence that the hard corrections to the inclusive top-quark pair-production at NNLO are indeed 
not large and it supports the estimate of our systematical uncertainty. 

To summarize, our approximate NNLO prediction leaves us with a rather small residual theo- 
retical uncertainty based on the scale variation. It is also worth stressing that the numerical impact 
of our theory improvements in Eqs. (fT2l) - (fT4l) and ([T8])-(l23l) is rather small, which again nicely 
illustrates the stability of the approximate NNLO predictions. In Tab. [D we compare with our pre- 
vious numbers [3,4] for the CTEQ6.6 [7] set at /% = 171 GeV and //^ = /i^ = m^. At Tevatron, we 
find no changes, as the cross section is entirely dominated by the ^^-channel at parton kinematics 
close to threshold. At LHC, there is a small net change of 2.5 pb in the prediction. Here, the effect 
of the improved NLO matching (in particular the In^jS term in Eq. (dH) depending on the constant 

in Eq. (fTTI) ) and the leading two-loop term fg^^^ in Eq. ([T3]) in the ^^-channel partially com- 
pensate. The numerical impact of the exact color decomposition at NLO (affecting the \n/3 term in 
Eq. ([141) ) is completely negligible. 
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LHC 


Tevatron 


o-Lo[pb] 


583.7 


5.820 


0-NL0[pb] 


877.4 


7.229 


C^NNLoLpb] (this work) 


923.0 


7.814 


C^NNL0[pb] (Refs. [3,4]) 


920.5 


7.810 



Table 1 : The LO, NLO and approximate NNLO prediction for the total cross section at LHC ( Vs = 14TeV) 
and Tevatron {^fs ^ 1.96TeV) using ^ 171 GeV, the PDF set CTEQ6.6 [7] and ^^=11^ = m,. For 
comparison we also give the previous numbers of Refs. [3,4]. 



As a central result of the current studies we quote our approximate NNLO prediction at LHC 
( VS = 14TeV) and Tevatron ( Vs = 1.96TeV) for a pole mass of = 173 GeV. For the MSTW 
2008 set [14] we have 



o-LHC = 887 pb 



^33 pb (scale) 
(TTev = 7.04 pb +2 pb (scale) 



pb(MSTW2008), 
IJIJpb (MSTW2008), 



and for CTEQ6.6 [7], 



CLHC 
(T'Tev 



874 pb !^3 



pb (scale) 



= 7.34 pb 



:J2Jpb (scale) 



pb(CTEQ6.6), 
pb (CTEQ6.6). 



(33) 
(34) 



(35) 
(36) 



Please note, that the IMSTW 2008 set [14] is based on a global analysis to NNLO in QCD while 
CTEQ6.6 [7] performs a fit to NLO only. Therefore, the two PDF sets return slightly different 
default values for the coupling constant or 5. While the LHC predictions of both sets are largely in 
agreement for these choices of a^, the difference in the Tevatron predictions can be attributed to 
differences in the parametrization of the light quark PDFs at large x. In addition, there is a system- 
atical uncertainty in Eqs. (l33])-(l36l) estimated to be (9(2%) due to unknown NNLO contributions. 



i.e. the exact expression for f'^^^^ in Eqs. (fT2l) - (fT4l) . One could of course argue that the accuracy 
of a given PDF set should match the accuracy of the theoretical prediction of the partonic cross 
section. However, we would like to disentangle the shift orginating from corrections to the hard 
parton scattering (which is the main subject of our paper) from PDF effects. Therefore we always 
use the same order in perturbation theory as far as the chosen PDFs are concerned (cf. Tab. [Hand 
also Appendix B of Ref. [3]). 

For applications, the mass dependence of the hadronic cross section ([T]) is conveniently param- 
eterized by the following simple fit formula 



cr(m,,/i) = a + bx + cx^ + dx^ + ex'^ + fx^ + gx^ , 



(37) 



with 1^=1^^= fij., X = (nit/GeV- 173) and the scale choices // = mt,2mt,mt/2. For the cross section 

at LHC (VS = 14TeV and ^fs = lOTeV) and Tevatron {^fs = 1.96TeV) all fit coefficients are 
listed in Tabs. [3] and H where we have used the PDF sets CTEQ6.6 [7] and IMSTW 2008 [14]. In 
the mass range 150GeV < < 220GeV, the accuracy of the fit is always better than 2.5%o. 

Let us briefly mention also other types of radiative corrections, which have not been considered 
here, e.g. in Eqs. (|33])-(IT7]). Within QCD these are bound state effects for the tt-pair near thresh- 
old [18, 19]. They affect the total cross section at LHC of the order (9(10) pb and, even more so. 
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differential distributions in the invariant mass of the top-quark pair. At Tevatron, due to the domi- 
nance of the ^^-channel in the color-octet configuration, they are negligible though. Precision anal- 
yses at the percent level naturally need to consider also the electro-weak radiative corrections at 
NLO [20-22]. Depending on the Higgs mass they cause a decrease relative to the LO cross section 
between 0(2%) for a light Higgs (w/, = 120GeV) and (9(2.5%) for a heavy Higgs (nih = lOOOGeV) 
at the LHC. This amounts to a negative contribution Acrgw - 0(lO- 15) pb. At the Tevatron, the 
electro-weak radiative corrections are almost zero for a light Higgs (m/, = 120GeV) and give a 
negative contribution of order (9(1%), i.e. Acrgw - (9(0.05) pb for a heavy Higgs (m/, = lOOOGeV). 



The top-quark mass in the MS scheme 

So far we have used the pole mass of the top-quark as a definition of the mass parameter. However, 
it is well-known that the concept of the pole mass has intrinsic theoretical limitations owing to 
the fact that the top-quark is a colored object. As such it does not appear as an asymptotic state 
of the 5 -matrix due to confinement. In other words the 5 -matrix does not have a pole in the 
top-quark channel. The impact of different mass renormalizations has been investigated in great 
detail in the context of top-quark mass measurements at a future linear collider where a precision 
of the order of a few hundred MeV is envisaged. In particular it has been shown that indeed the 
conceptual limitations of the pole mass lead to a poorly behaved perturbative series. A class of 
alternative mass definitions, so-called short distance masses, offer a solution to this problem, e.g. 
the IS -mass or the potential subtracted (PS) mass (see e.g. Ref. [41]). 

In the following we study the impact of the conversion from the pole mass scheme to the MS 
scheme (see Ref. [23-25] and references therein) for the total cross section of top-quark hadro- 
production. This is a novel feature and, in principle, the cross section in terms of the MS mass can 
be used for a direct measurement of the running mass at a high scale. This is similar to the case 
of Z?-quark production at LEP (see Refs. [42,43] and [44-47]). Let us first describe briefly how 
we translate the predictions for the total cross section from the pole mass to the MS mass scheme. 
The starting point is the well-known relation between the pole mass and the MS mass m(//^) to 
NNLO: 

nit = m(jj.^)(^l+as(Mr)d^^^ +as(iUrf'd^^^^ , (38) 

with Us = a's^ jn (i.e. five active flavors) and coefficients d^^\ which in general depend on the 
ratio fj.^^/m(Hj.)^, 

d^'^ = (39) 
,r2> 307 ^ 2,^ 1 509, 47, 2 

- + ^^2 + ^^mo-,) + + ^ Z ^("^'/m,) . (40) 

Here nf denotes the number of light flavors and = ln(//^^/m(jU^)^). The function A(mi/mt) 

accounts for all massive quarks m, lighter than the top-quark. For all light quarks we set m, = so 
the sum in Eq. (|40|) vanishes. Note also that the decoupling of the top-quark is assumed to be done 
at the scale of the MS mass in(/u^). 

Let us start by making the mass dependence in the total cross section manifest order by order 
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in perturbation theory. For the pole mass we have through NNLO 



(41) 



! = 



Next, we use the relation (l38l) above to convert from the pole mass to the MS mass m(m). For 
simplicity we abbreviate m = m(m) and obtain 



(T : 



( 



1=0 



1=1 



ii-l) 



(m) 



(42) 



We note that the coefficients S'"* have to be evaluated for /i,. = m (corresponding to the scale of 
as). Thus, the task in Eq. (I42l) amounts to determine the derivatives of the cross sections cr^'^ with 
respect to the mass. To do so in practice we have chosen the following approach. For all coefficients 
(T^'^ we use the ansatz of Eq. (l37l) to parametrize the mass dependence. More precisely we evaluate 
the hadronic cross section order by order in perturbation theory for a fixed renormalization and 
factorization scale. Then, varying the top-quark mass (in the pole mass scheme) and performing 
a fit similar to what has been discussed before in Eq. (|37l) we obtain the total cross section in the 
following form: 



N 



(43) 



;=0 



k=0 



where c^'^ denote the (order dependent) fit coefficients. is the order of the polynomial in trit 
{N = 6'm Eq. (ITT]) ) and mo is our fixed reference mass (taken to be 173 GeV in Eq. (|37l) ). Since all 
dependence on the pole mass is manifest, it is now a straightforward exercise to convert to the 
MS mass m and to perform the derivatives in Eq. (|42)) . 



A' 



cr 



(=0 



k=0 

1 



.(i-l) 



(44) 



/=i 



+6i^2^k{k-\)m^{m-mo)'^ ^ 



niQ—m 



JO) 

-1 



Tr:2 



JO) I 



If the expansion point mo is chosen to be the MS mass m, Eq. (|44|) simplifies considerably and 
the truncation of the power series in m to first (second) order is exact at NLO (NNLO). Generally 
though, for applications, it is of some advantage to keep mo at a fixed numerical value and to 
rely on the fact, that our ansatz (|43l) with a polynomial of high enough degree approximates all 
coefficients cr^'^ and their first two derivatives sufficiently well. As discussed below Eq. (ITTl) . the 
choice N = 6 achieves per mille accuracy in the phenomenologically interesting range. We have 
also checked that the choices mo = m and niQi^m yield the same result. 

We stress again, that we have fixed /i^ = m in Eq. (I44l) . However, it is also possible to restore the 
complete renormalization scale dependence using the well-known relation for the running coupling 



a,(m) = a,(/z^)(l +4n^as(M^)L^j3o + (4n^fas(M^f(j3iL^ +^0^1)) ' 



(45) 
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Figure 4: The scale dependence of the total cross section with the top-quark mass in the MS scheme at 
m = 163GeV at LO (red), NLO (green) and approximate NNLO (blue). The dashed lines denote the ji^ -m 
for the factorization scale, the solid lines the maximal deviations for //^ e [m/2,2m] and /i^ = m/2,m and 2m. 

We use the MSTW 2008 PDF set [14] at LHC with - 14TeV (left) and Tevatron with y/S - 1.96TeV 
(right). The vertical bars indicate the size of the scale variation in the standard range [m/2,2rn]. 

with = ln(/i^^/m^) and fio and fii given in Eq. (I24l) . To summarize, Eq. (I44l) represents an 
explicit expression for the total cross section of top-quark hadro-production with the top-quark 
mass defined in the MS scheme. 

Let us illustrate the phenomenological consequences of the MS mass for predictions at Tevatron 
and LHC. In Fig. |4] we plot the scale dependence of the total cross section again at the various 
orders in perturbation theory. The value of m = 163GeV roughly corresponds to a pole mass of 
nif = 173 GeV and we choose three (fixed) values for the factorization scale fi^ = m/2,m and 2m. 
The band to the left denotes the maximum and the minimum values of ju^ 6 [m/2, 2m] for the three 
choices of fi^ according to Eq. (|3T1) . cf. the contour plot in Fig. [T] for the pole mass. We observe a 
great stability with respect to scale variations when including higher order perturbative corrections 
through NNLO. Remarkably, at Tevatron, the scale variation at NNLO is even reduced further by 
more than a factor of two compared to the result in the pole mass scheme. 

Next, in Fig. [5] we show the mass dependence of the total cross section employing the MS 
mass definition and performing the same scale variation as above, i.e. /i^ = m/2,m and 2m and 
yu^ 6 [m/2, 2m]. Upon adding the higher order perturbative corrections we observe as a striking 
feature the extremely small numerical effect of the radiative corrections. E.g. for m = 163 GeV at 
Tevatron, we find the effect of the NLO corrections to be only 1.5% and even less (0.9%) for the 
approximate NNLO results. Also for the LHC, we observe a much faster convergence of the per- 
turbative expansion when using the MS mass. The NLO (approximate NNLO) corrections amount 
to 31% (4%) at m = 163 GeV which is roughly half of the size of the corrections in the pole mass 
scheme. This demonstrates an extremely good stability of the perturbative series in the MS mass 
scheme. We can understand this behavior qualitatively by looking at the mass dependence of the 
scaling functions in Eqs. (|4l)-(l9l). We find e.g. dmff^^ - (1 -/3^)//3, which implies sizably enhanced 
contributions near partonic threshold, i.e. precisely in the region which contributes dominantly 
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Figure 5: The mass dependence of the total cross section for the MS mass m at NLO (green) and approxi- 
mate NNLO (red) order with the scale variation in the range ju^ € [m/2,2m] and ju^ = m/2,m and 2rn for the 

MSTW 2008 PDF set [14] at LHC with y/S - 14TeV (left) and Tevatron with yfs = 1.96TeV (right). The 
value for the Tevatron cross section is taken from Ref. [8]. 



in the convolution with the parton luminosities, cf. Eq. ([T]). This observation is yet another ar- 
gument in favor of the phenomenological importance of our approximate NNLO predictions in 
Eqs. (Hl-dlll). 





m [GeV] 


mt [GeV] 


LO 

NLO 

NNLO 


159.2^3.5 

159.8!3;3 
160.0+3-3 


159.2+3-5 
165.8:3J 



Table 2: The LO, NLO and approximate NNLO results for the top-quark mass in the MS scheme (m) and 
the pole mass scheme (nit) for the measured cross section of cr = 8.18pb at Tevatron [8]. The uncertainties 
in the Table reflect the quoted experimental uncertainties. 

A different way to address the issue of perturbative stability is the extraction of the IMS mass 
from the total cross section as measured at Tevatron. Ref. [8] quotes a value with a combined 
uncertainty of cr = 8. 1 8+q pb for a top-quark mass m/ = 170 GeV along with a (weak) dependence 
on the value of the mass, e.g. cr = 7.99 pb for the latest world average [13], = 173.1+| 2GeV. 
Using the measured value of cr = 8.18pb which is consistent with the theory predictions of Refs. [3, 
4] (and with this work, see Tab. O we extract the MS mass m order by order. As mentioned 
earlier we use the same NNLO PDF set of IMSTW 2008 [14] independent the order of perturbation 
theory and the results at LO, NLO and approximate NNLO are given in Tab. l2l The value of 
m = 160.0+3-2 GeV represents to the best of our knowledge the first direct determination of the 
running top-quark mass from experimental data. For comparison, we also quote the values of the 
pole mass at the respective order extracted in the same way by directly comparing the theory 
prediction with the measured cross section. Alternatively, we can also convert the MS mass value 
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back to the pole mass scheme with the help of Eq. (|38l ). Our NNLO value for m corresponds to 
mt = 168.9]J;2 4GeV, which constitutes a theoretically well-defined determination of the pole mass 
and is also (within the experimental uncertainties) in agreement with the world average [13] of 
nit = IVB.l^l^j gGeV. To summarize, the MS mass scheme is distinguished by the great stability 
in the value of the extracted top-quark mass. This feature has been studied in the past in detail 
for processes at a future linear collider [41] and our observation is also in agreement with recent 
considerations based on the renormalization group flow for heavy quark masses [48]. 

Summary 

In this Letter, we update and extend the predictions of Refs. [3,4] for the cross section of top-quark 
hadro-production at LHC and Tevatron. We have applied some improvements in the threshold 
approximation for the two-loop scaling functions ([T2l)-([T4l) as described in the text. We provide 
new and precise parametrizations in Eqs. ([T])-® and ([T8l)-(|23l) for all scaling functions that can 
be determined exactly. All fit functions are documented in the Appendix. Moreover, we have 
performed the independent variation of the renormalization and the factorization scale with the 
help of Eqs. (|29|) . (l30l) . As a novel aspect, in addition to the conventionally used pole mass we 
provide predictions for the total cross section employing the MS definition for the mass parameter. 
The central result is Eq. (I44l) . 

Our main phenomenological results are the parametrizations in Eq. (l37l) . Tabs. [3] and |4] and the 
cross sections in Eqs. (l33])-(l36l) for the pole mass. The diff'erences with respect to our previous 
numbers are quite small, though, see Tab. [TJ The theory uncertainty according to Eq. (|3TI) defined 
by exploring the (/i^,/i^)-plane in the standard range yu^,//^ 6 [m;/2,2mf] does not differ signifi- 
cantly from the case of fixed scales ji^ = jj.^. We have also addressed the residual systematical 
uncertainty due to the threshold approximation and we have quantified the effect of other higher 
order corrections, such as electro-weak or QCD bound state effects. The most interesting aspect 
of our phenomenological studies consists of the conversion to the MS mass scheme in Figs. |4] and 
[51 The cross section predictions with the MS mass definition exhibit a greatly improved pattern of 
apparent convergence for the perturbative expansion and very good stability with respect to scale 
variations. This leads to very stable values for the extracted mass parameter m as given in Tab. |2l 
In particular we find 

m= leO.O^^-^GeV. (46) 

This is the first direct determination of the running top-quark mass from top-quark pair-production. 
The corresponding value for the pole mass derived from Eq. (|46|) reads 

= 168.9!^jGeV, (47) 

which is consistent with current world average [13], rut = n3.l^^^^l^ GeV. Altogether, this provides 
substantial support in view of the reliability of our approximate NNLO numbers. We believe 
that the QCD radiative correction for top-quark pair-production at hadron colliders are well under 
control. 
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where all threshold logarithms ln(/3) and the Coulomb corrections (~ 1 /J3) are exact. The fit func- 
tions are given in Eqs. (|A.7| )- (|A.9I) and all parameters of the fit are listed in Tabs. [5]-(7] The fits to 
the scaling functions fij^\ fjj^^ in Eqs. (IA.1I) - (IA.6I) are, in general, accurate at the per mille level. 
Exceptions are regions close to zero, which is not surprising. There we retain an accuracy better 
than one percent. 

Fortran subroutines with the parametrizations of all scaling functions and the coefficient in 
Tabs. [5H7] are available from the authors upon request. 

h(fi,ai,...,an) = a{/3^ + a2/3^ +a3/3* + a4/3^ (A.7) 
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+a9yS^ In^yS + a 1 oyS^ In^ /3 + a 1 1/3 Inp + a 1 2/3 In^ p + a 1 3jS^ Inp 
+ai4/3^ In^p + a 15/3-^ Inp + ai^p^ In^p + ai']0^ Inp , 

+a4ye'^ liiyS + asp^ \nl3 + a^fi^ ln/3 
+a-]fi^p Inp + asySV In^p + ag/^V Inp 
+a 1 oyS^p In^ p + a 1 i0^p In p 

+a\2^p In^p + ai3/3^p In^ p + a^jS^p In^p + aisfi^p In^p , 
hf^(J3,ai,...,ai8) = aip^ + a20^ + a^p^ + a^p^ + asp' (A.9) 

+a6>e'^ InyS + a7/3^ InjS + asyS^ ln/3 + 09/3^ ln/3 
+aio/3-^ Inp + aii/3^ In^p + a 12/3"^ Inp + ai30^\n^ p 
+auP^ Inp + ai5)S^ In^p + a i6/3^ Inp + a nyS^ In^p + a i8/3^ Inp . 
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-4.45624x10" 


-3 


4.81468x10" 


-5 


-4.39515x 10" 


-7 


2.28915x10" 


-9 


crip. = m,) + Acr 


9.02378 X 10^ 


-2.40942 X 10' 


3.81862 X 10"' 


-4.71806x10" 


-3 


5.10568x10" 


-5 


-4.67281 X 10" 


-7 


2.43820x10" 


-9 


cr{p = m,) - Acr 


8.46479 xlO^ 


-2.29441 X 10' 


3.66298 X 10"' 


-4.53009x10 


-3 


4.88109x10" 


-5 


-4.43754 X 10" 


-7 


2.30362x10" 


-9 



LHC - lOTeV, CTEQ6.6 



crip 


= mt) 


3.96877 X 10^ 


-1.12077 X 10' 


1.85352 X 10"' 


-2.36659x10" 


-3 


2.62800x10" 


-5 


-2.44841 X 10" 


-7 


1.28959x10" 


~9 


crip 


= m,/2) 


3.97 124 xlO^ 


-1.11889X 10' 


1.84706x10"' 


-2.35501 X 10" 


-3 


2.61183x10" 


-5 


-2.42989 X 10" 


-7 


1.27805x10" 


-9 


crip 


= 2m,) 


3.79852 xlO^ 


-1.07358x10' 


1.77667x10"' 


-2.26977 X 10" 


-3 


2.52223x10" 


-5 


-2.35220 X 10" 


-7 


1.24016x10" 


-9 


crip 


= m,) + Acr 


4.15125 xlO^ 


-1.15947 X 10' 


1.90285x10"' 


-2.41772x10" 


-3 


2.67843x10" 


-5 


-2.49480 X 10" 


-7 


1.31488x10" 


-9 


a-(p 


= m,) - Acr 


3.78628 xlO^ 


-1.08207 X 10' 


1.80416 X 10"' 


-2.31532x10" 


-3 


2.57769x10" 


-5 


-2.40323 X 10" 


-7 


1.26569x10" 


-9 



Tevatron - 1.96TeV, CTEQ6.6 



aip 


= mt) 


7.34317x10" 


-2.27486x10"' 


3.94086x10" 


-3 


-5.22302x10" 


-5 


6.09497 X 10 


-7 


-5.99414 X 10" 


-9 


3.27925 X 10- 


ll 


crip 


= m,/2) 


7.58312x10" 


-2.34571x10"' 


4.05822x10" 


-3 


-5.37018x10" 


-5 


6.25408 X 10" 


-7 


-6.13901 X 10" 


-9 


3.35467 X 10" 


11 


crip 


= 2m,) 


6.96303 X 10" 


-2.15748x10"' 


3.73128x10" 


-3 


-4.93012x10" 


-5 


5.73218x10" 


-7 


-5.62092 X 10" 


-9 


3.07038 X 10" 


11 


crip 


= m,) + Act 


7.75854x10" 


-2.42254x10"' 


4.23665 X 10" 


-3 


-5.65955x10" 


-5 


6.63296x10" 


-7 


-6.52935x 10" 


-9 


3.57062 X 10" 


11 


crip 


- nit) - Act 


6.92780x10" 


-2.12718 X 10"' 


3.64506x 10" 


-3 


-4.78628x10" 


-5 


5.55679x10" 


-7 


-5.46023 X 10" 


-9 


2.99003 X lo- 


11 



Table 3: Fit coefficients to Eq. (l37l) for crip = mt,2mt,m,l2) and crip = m,) ± Act for the PDF set CTEQ6.6 [7] and the colliders LHC and Tevatron. 



a[ph] b[ph] c[pb] d[ph] e[ph] f[ph] g[ph] 



LHC ^ 14TeV, MSTW 2008 NNLO 



o-Qj. 


= mt) 


8.87496x102 


-2.38344 X 10' 


3.78224 xlO^' 


-4.66307 X 10" 


-3 


5.02155x10" 


-5 


-4.56910 X 10" 


-7 


2.37374x10" 


-9 


a(ji 


= m,/2) 


8.85530x102 


-2.37387 X 10' 


3.76203x10-' 


-4.63331x10" 


-3 


4.98411x10" 


-5 


-4.52980 X 10" 


-7 


2.35137x10" 


~9 


o-Qj. 


= 2m,) 


8.54052 xlO^ 


-2.29566 X 10' 


3.64547 X 10-' 


-4.49661x10" 


-3 


4.84539x10" 


-5 


-4.41574 X 10" 


-7 


2.29907 X 10" 


-9 


crip. 


= nit) + Ao" 


9.02902 xlO^ 


-2.41907 X 10' 


3.83190x10"' 


-4.71808x10" 


-3 


5.07642x10" 


-5 


-4.61804 X 10" 


-7 


2.39989x10" 


-9 




= nit) - ^cr 


8.72090x10^ 


-2.34783 X 10' 


3.73257x10-' 


-4.60776x10 


-3 


4.96661x10 


-5 


-4.52297 X 10" 


-7 


2.35168x10" 


-9 



LHC - lOTeV, MSTW 2008 NNLO 



crip 


= nit) 


4.03219 xl02 


-1.13904X 10' 


1.88177 X 10-' 


-2.39835x10" 


-3 


2.65811x10" 


-5 


-2.47337 X 10" 


-7 


1.30217 X 10" 


^9 


crip 


= m,/2) 


4.03439 xlO^ 


-1.13695x 10' 


1.87488x10"' 


-2.38625x10" 


-3 


2.64067 X 10" 


-5 


-2.45191 X 10" 


-7 


1.28831x10" 


-9 


crip. 


= Intf) 


3.86012 xlO^ 


-1.09154x10' 


1.80486x10-' 


-2.30194x10" 


-3 


2.55275x10" 


-5 


-2.37710 X 10" 


-7 


1.25272x10" 


'9 


crip. 


= nit) + Ao" 


4.11912x10^ 


-1.16047 X 10' 


1.91287x10-' 


-2.43344x10" 


-3 


2.69297 X 10" 


-5 


-2.50369x 10" 


-7 


1.31793x10" 


-9 


crip 


= Mt) - Acr 


3.94526 xlO^ 


-1.11761 X 10' 


1.85066x 10-' 


-2.36310x10" 


-3 


2.62315x10" 


-5 


-2.44419 X 10" 


-7 


1.28819 X 10" 


-9 



Tevatron - 1.96TeV, MSTW 2008 NNLO 



crip 


= nit) 


7.04217x10" 


-2.18800X 10-' 


3.80366x10" 


-3 


-5.06795 X 10" 


-5 


5.96308x10 


-7 


-5.92150x 10" 


-9 


3.26369x10" 


ll 


crip 


= mt/2) 


7.27746x10" 


-2.25794x10-' 


3.92060x10" 


^3 


-5.21816x10" 


-5 


6.12661x10" 


-7 


-6.05762 X 10" 


-9 


3.32360x 10" 


11 


crip 


= Inif) 


6.67970x10" 


-2.07517x10-' 


3.60070x10" 


-3 


-4.78294x10" 


-5 


5.60141x10" 


-7 


-5.52680x 10" 


-9 


3.02883x 10" 


11 


crip 


= nit) + Act 


7.18407x10" 


-2.22429x10-' 


3.85513x 10" 


-3 


-5.13171x10" 


^5 


6.05482x10" 


-7 


-6.04397 X 10" 


-9 


3.34549x10" 


11 


crip 


- nit) - Act 


6.90028 X 10" 


-2.15171 X 10-' 


3.75214 X 10" 


-3 


-5.00394x10 


-5 


5.87170x10" 


-7 


-5.80214 X 10" 


-9 


3.18553 X lo- 


11 



Table 4: Same as Tab.[3]for the PDF set MSTW2008 [14] at NNLO. The PDF uncertainty Act has been obtained with the 68% confidence level set. 



i 


Jqq 


^(21) 

Jqq 

hi Ci 


.(22) 
Jqq 

hi Ci 


1 

1 


0.07120603 


-0.15388765 


-0.07960658 


0.37947056 


-0.00224114 


2 


-1.27169999 


4 85226571 


50111294 


-4 25138041 


02685576 


3 


1.24099536 


-7.06602840 


-0.09496432 


2.91716094 


-0.01777126 


4 


-0.04050443 


2.36935255 


-0.32590203 


0.94994470 


-0.00626121 


5 


0.02053737 


-0.03634651 


-0.02229012 


0.10537529 


-0.00062062 


6 


-0.31763337 


1.25860837 


0.23397666 


-1.69689874 


0.00980999 


7 


-0.71439686 


2.75441901 


0.30223487 


-2.60977181 


0.01631175 


8 


0.01170002 


-1.26571709 


0.13113818 


-0.27215567 


0.00182500 


9 


0.00148918 


-0.00230536 


-0.00162603 


0.00787855 


-0.00004627 


10 


-0.14451497 


0.15633927 


0.08378465 


-0.47933827 


0.00286176 


11 


-0.13906364 


1.79535231 


-0.09147804 


-0.18217132 


0.00111459 


12 


0.01076756 


0.36960437 


-0.01581518 


-0.04067972 


0.00017425 


13 


0.49397845 


-5.45794874 


0.26834309 


0.54147194 


-0.00359593 


14 


-0.00567381 


-0.76651636 


0.03251642 


0.08404406 


-0.00035339 


15 


-0.53741901 


5.35350436 


-0.25679483 


-0.51918414 


0.00363300 


16 


-0.00509378 


0.39690927 


-0.01670122 


-0.04336452 


0.00017915 


17 


0.18250366 


-1.68935685 


0.07993054 


0.15957988 


-0.00115164 



Table 5: Coefficients for fits of the qq scaling functions. 
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.(10) 


^(21) 




.(22) 

Jgq 


/ 


Ui 


hi 


Ci 


h- 




1 


-0.26103970 


-0.00120532 


0.00003257 


-0.00022247 


0.00001789 


2 


0.30192672 


-0.04906353 


0.00014276 


0.00050422 


0.00000071 


3 


-0.01505487 


-0.20885725 - 


-0.00402017 


-0.02945504 


-0.00020581 


4 


-0.00142150 


-13.73137224 


0.06329831 


0.34340412 


0.00108759 


5 


-0.04660699 


14.01818840 - 


-0.05952825 


-0.31894917 


-0.00086284 


6 


-0.15089038 


-0.00930488 


0.00002694 


0.00009213 


0.00000010 


7 


-0.25397761 


-0.52223668 


0.00159804 


0.00690402 


0.00001638 


8 


-0.00999129 


-4.68440515 


0.01522672 


0.07847233 


0.00022730 


9 


0.39878717 


-7.61046166 


0.02869438 


0.16042051 


0.00045698 


10 




1.36687743 - 


-0.00875589 


-0.05186974 


-0.00025620 


11 


-0.14178346 


1.84698291 - 


-0.00800271 


-0.03861021 


-0.00016026 


12 


0.01867287 


-7.26265988 


0.04043479 


0.21650362 


0.00070713 


13 


0.00238656 


-4.89364026 


0.01965878 


0.10137656 


0.00034937 


14 


-0.00003399 


11.04566784 - 


-0.05262293 


-0.28056264 


-0.00072547 


15 


-0.00000089 


4.13660190 - 


-0.01457395 


-0.08090469 


-0.00025525 


16 


0.00000000 


-6.33477051 


0.02314616 


0.13077889 


0.00034015 


17 


0.00000000 


-1.08995440 


0.00291792 


0.01813862 


0.00006613 


18 


0.00000000 


1.19010561 - 


-0.00220115 


-0.01585757 


-0.00006562 



Table 6: Coefficients for fits of the gq scaling functions. 
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/ 


a, 


Jgg 

b, 




All) 
Jgg 

bi Ci 


1 

1 


-8.92563222 


-4.18931464 


0.12306772 


0.01222783 


-0.00380386 


2 


149.90572830 


82.35066406 - 


-2.75808806 


-0 77856184 


08757766 


3 


-140.55601420 


-87.87311969 


3.19739272 


1.33955698 


-0.10742267 


4 


-0.34115615 


9.80259328 - 


-0.56233045 


-0.59108409 


0.02382706 


5 


-2.41049833 


-1.12268550 


0.03240048 


0.00248333 


-0.00099760 


6 


54.73381889 


29.51830225 - 


-0.92541788 


-0.23827213 


0.02932941 


7 


90.91548015 


48.36110694 - 


-1.57154712 


-0.38868910 


0.04906147 


8 


-4.88401008 


-7.06261770 


0.35109760 


0.28342153 


-0.01373734 


9 


-0.17466779 


-0.08025226 


0.00227936 


0.00010876 


-0.00006986 


10 


13.47033628 


7.01493779 - 


-0.21030153 


-0.03383862 


0.00658371 


11 


22.66482710 


15.00588140 - 


-0.63688407 


-0.29071016 


0.02089321 


12 


4.60726682 


3.84142441 - 


-0.12959776 


-0.11473654 


0.00495414 


13 


-67.62342328 


-47.02161789 


1.91690216 


0.98929369 


-0.06553459 


14 


-9.70391427 


-8.05583379 


0.26755747 


0.24899069 


-0.01046635 


15 


65.08050888 


47.02740535 - 


-1.86154423 


-1.06096321 


0.06559130 


16 


5.09663260 


4.21438052 - 


-0.13795865 


-0.13425338 


0.00551218 


17 




-14.99599732 


0.58155056 


0.35935660 


-0.02095059 



Table 7: Coefficients for fits of the gg scaling functions. 
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